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ABSTRACT 



^ ■ A model of a long optical communication line consisting of alternating segments with 



anomalous and normal dispersion, whose lengths are picked up randomly from a certain 
interval, is considered. As the first stage of the analysis, we calculate small changes 
of parameters of a quasi- Gaussian pulse passing a double-segment cell by means of the 
variational approximation (VA) and approximate the evolution of the pulse passing many 
cells by smoothed ODEs with random coefficients, which are solved numerically. Next, we 
perform systematic direct simulations of the model. Results are presented as dependences 
of the pulse's mean width, and standard deviation of the width from its mean value, on 
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the propagation distance. Averaging over 200 different realizations of the random length 
set reveals slow long-scale dynamics of the pulse, frequently in the form of long-period 
oscillations of its width. It is thus found that the soliton is most stable in the case of 
the zero path-average dispersion (PAD), less stable in the case of anomalous PAD, and 
least stable in the case of normal PAD. The soliton's stability also strongly depends on 
its energy, the soliton with small energy being much more robust than its large-energy 
counterpart. 



1 Introduction 

The use of dispersion management (DM) for improvement of pulse transmission in long 



optical- fiber links has attracted a great deal of attention, see, e.g., Refs. []I|-|1T2|| and 
references therein. Almost all the works on this topic published thus far were dealing 
with conditions for stationary propagation of a pulse in a link consisting of periodically 
alternating fixed-length pieces of the fiber with anomalous and normal dispersion. How- 
ever, really exiting terrestrial optical communication webs are patchwork systems, which 
include links with different values of the dispersion and, what is especially important, very 
different lengths. If the distribution of the lengths may be assumed random, upgrading 
the patchwork systems by means of DM makes it necessary to consider transmission of 
pulses in long communication lines subject to random DM. Besides the obvious signifi- 
cance for the applications, this issue if also of considerable interest by itself, in the context 
of soliton propagation in strongly inhomogeneous media (see, e.g., the book [|TB1). The 



objective of this work is to develop both semi-analytical variational and fully numerical 
approaches to the description of pulses' dynamics in a random-DM system. The varia- 
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tional approximation (VA), which was first applied to the DM models in |14 |, has then 
become a commonly accepted tool for theoretical consideration of DM solitons, see, e.g., 
Refs. |, |, |]. 

Most promising results, viz., solitons that are most robust against random variation 
of the fiber-segment lengths, will be found in this work for the case of zero path-average 
dispersion (PAD) /3q. Note that the case of zero or slightly normal (/3q > 0) PAD has 
recently attracted special attention in the context of regular (periodic) DM, after it had 
been demonstrated that stable propagation of pulses is possible in this case too j|, || § . 

There are two different natural models of random DM. In the two-lengths (2L) model, 
the lengths of the alternating anomalous- and normal-dispersion pieces are selected ran- 
domly from a certain interval independently from each other. In the other, one-length (1L) 
model, equal lengths of the anomalous- and normal- dispersion pieces inside each DM cell 
are selected randomly from the same interval. In the 2L model, fio is the global average 
value, while in the 1L model (3q is the mean dispersion in each cell. Final results are quite 
similar for both models, therefore in this work they will be presented for the 1L version. 

Evolution of the soliton is a random setting is quantified by a dependence of its width 
on the propagation distance. In the random-DM links, the soliton's width demonstrates 
rapid erratic oscillations. However, averaging over a large number of the random-DM 
realizations, and taking average deviations from the mean values, we obtain smooth re- 
sults which can be understood. Roughly, the main conclusions are that robustness of a 
randomly dispersion- managed soliton crucially depends on the value of PAD and on the 
soliton's energy. First, the soliton's degradation (decay) is fast at (3q > 0, slower at [3q < 
(recall this corresponds to the normal and anomalous dispersion, respectively), and much 
slower at /3 = 0. Second, degradation is much slower for solitons with small energy than 



for those with large energy. Both these findings can be understood if one notices that, in 
the exactly solvable linear model (corresponding to the small-energy limit) with random 
DM, a chirped Gaussian pulse is chaotically vibrating with no systematic degradation 
(spreading out), provided that PAD is exactly zero (see, e.g., Ref. f|). 

A preliminary version of the present work was published in ||. Recently, another 
random-DM model was also considered by Abdullaev and Baizakov || (see also [|l"5|j ). 
However, in that work emphasis was made on a different model, in which the local values 
of the dispersion, rather than the fiber-segment lengths, were subject to a random dis- 
tribution (the case of the randomly distributed lengths was only briefly touched within 
the framework of VA in Ref. |J). Unlike the results reported in the present work, sharp 
differences between robust and degrading soliton propagation regimes were not found in 
the model considered in Ref. ||, 

The rest of the paper is organized as follows. In section 2, the random-DM model is 
formulated. In section 3, we recapitulate basic results of VA for the periodic-DM model, 
and then develop VA for the random-DM one; the approximation is based on the Gaussian 
ansatz for the pulse. In section 4, results of systematic numerical simulations of the VA- 
generated equations with randomly varying coefficients are displayed. Section 5 presents 
results of direct numerical simulations of the underlying nonlinear Schrodinger (NLS) 
equation with random DM, which are, generally, found to be in accord with predictions 
of VA. Section 6 concludes the paper. 
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2 The random-dispersion-management model 

We start with a model which is a straightforward generalization of the well-known DM 
schemes based on the NLS equation 

iu z - ^P(z)u TT + |-u| 2 m = 0, (1) 



which is written in the standard notation ||17|1 . The dispersion coefficient is taken in the 



form (3{z) = (3q + /3i(z), where f3 Q is its average value (PAD), and (3\{z) is a variable part 
with the zero average, which is taken in the following form inside the n-th DM cell: 

{/?_ , z n < z < z n + , 
(2) 
P+ , z n + LH < z< z n + LH + LW = Zn+1 . 

Here, /3_ +/3q < and (3+ + /3 > are, respectively, the dispersion coefficient in the 

anomalous- and normal- dispersion fibers and L^' are lengths of the corresponding pieces. 

In the case of periodic DM, the lengths are the same in all the cells, while in the 

case of random DM, they vary stochastically from a cell to a cell. The condition that 



the average value of Pi(z) is zero implies that /3_L( _ ) + P+L^ = 0, the overbar standing 
for the averaging. We assume that random values of both and are distributed 
uniformly within a certain interval, so that their mean values are equal. Then, the above 
condition requires that = (3 +) which will be assumed to hold. 
The model ([!]) conserves the total energy E, which we define as 

/ r+oo 

E=J2/ir \u{r)\ 2 dT. (3) 

i-oo 



Following many earlier works P]]- HTDfl , the model does not include losses and gain, pre- 
suming that they are mutually compensated at a smaller scale. The analysis presented 
below can be extended to include the losses and gain, but this is postponed to another 



work, as it is necessary first of all to understand properties of the conservative model. 
Other effects which are not included into the present work but may be relevant in certain 
cases, as it is known that they may help to stabilize the transmission of pulses in the long 
periodic-DM link, are the third-order dispersion [H| and filtering |TT |. 



We are concerned with the case when DM is not too weak, hence the shape of the 
pulse is well-known to be close to a Gaussian flfl-|10]|. Accordingly, VA may be based 
on the Gaussian ansatz [0- ||. We will here follow a version of the Gaussian-based VA 
developed (for periodic DM) in detail in ||. Using the scaling invariance of Eq. ([!]), the 
following normalizations are adopted in Ref . || , 

L i-) + L (+) = i ; \p_\ £(-) = p +L (+) = i. (4) 
In the present work, we apply these normalizations to the mean values of the random 



lengths. Since our model assumes = L(+), Eq. (H) yields = 1/2, and \/3±\ = 2. 
To comply with the former normalization, we choose the interval from which the random 
lengths are picked up as 

0.1<L<0.9. (5) 

The minimum length 0.1 is introduced here because, in reality, the length can be neither 
very large (say, larger than 200 km) nor very small (shorter than 20 km). 



3 The variational approximation 

3.1 The Gaussian ansatz 

The strong-DM regime implies that, locally, the dispersion is much stronger than the 
nonlinearity, and that \/3q\ <C |/3=p|. In the zero-order approximation, completely neglecting 
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the nonlinearity, one has an exact Gaussian solution to the linearized equation (|]) ||, 

(6) 



u = = exp 

It 2 + 2iA(z) 



t 2 



+ ^c 



t 2 + 2i/\(z) 

Here, Pq and r are, respectively, the peak power and minimum width of the pulse, — A(z) 
= — Ao + H n f3(z')dz' is the accumulated dispersion defined inside the n-th DM cell, and 
Ao and <fi are real constants. 

The parameter Tq 2 in the expression (^) is proportional, with regard to the normaliza- 
tions (HD, to the well-known DM strength, which (in the case of more general normaliza- 
tions) is defined as @ S = (|/3_| LH + (3 + L^) /r| WHM , where t F whm is the full width at 



half maximum (FWHM) of the pulse |TJ]] at the midpoint of the anomalous-dispersion seg- 
ment, where the pulse is narrowest. Applying the notation adopted here, we replace L^' 
in the definition of S by the above mean values 1/2, and also insert \/3±\ = 2, obtaining 

5 = 1.443/r 2 . (7) 

The strength S is the most important characteristic of the DM schemes. It determines 
their basic properties, which virtually do not depend on other parameters (such as, e.g., 
L + /L_), provided that S is fixed. In particular, detailed numerical simulations reveal 
that stable DM pulses do not exist at S > S max ~ 10, the propagation at zero or slightly 
normal PAD is possible is S > S CT ~ 4 j|, |5|, || , and strongest suppression of the interaction 



between the pulses is attained at S ~ 1.6 [12 



The exact solution (|6|) for the linear model will be used below as an ansatz on which 
VA for the nonlinear model is based. In most other versions of VA [RJ, the Gaussian 
ansatz is also used, but in a different form, 

u = a(z) exp \-Tq/W 2 (z) + ib(z)r 2 +i<f>] . (8) 



The complex amplitude a(z) and real width W(z) and chirp b(z) introduced in this 
expression are related to parameters of the ansatz (H) as follows: 



-i 



a 2 {z) = r 2 P [r 2 + 2tA{z)\ , W{z) = T^y/ifi + 4A 2 (z) , b(z) = 2A [r 4 + 4A 2 (^ 

(9) 

Using the VA technique, one can derive equations for the nonlinearity-induced evo- 
lution of the parameters P, r , and A , that were constant within the framework of the 
exact Gaussian solution in the absence of the nonlinearity ||. First, in accord with the 
energy conservation, we obtain P To = E — const (this coincides with the conserved 
energy defined by Eq. (||) ), and then 

dr _ V2Et A(z) rfAo E [AA 2 (z) — Tq] 

dz W\z) ' dz ~ P0+ 2V2W 3 (z) ' 1 J 

Because the average dispersion is small in the DM regime, it is also treated (to derive the 
second equation in (|10|)) weak perturbation. 

The changes of the parameters r and A per one DM cell can be calculated as 

where § stands for the integration over a full cell, from z = z n to z = z n + L^> + 
In the spirit of the perturbation theory, the changes (|TTJ) are assumed small, hence, 
calculating the integrals in Eq. (|TTf) , it is sufficient to take into regard only the rapid 
variation of A(z), while tq and Ao are treated as constants. 

3.2 Revisiting the case of the periodic dispersion management 

To develop VA for the random-DM system, it is first necessary to recapitulate basic results 
for the usual periodic case. In that case, to obtain conditions providing for the stationary 
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transmission of the Gaussian pulse in the long DM line, one equates to zero 5tq and 5Ao, 
evaluated as per Eq. (|TT|), with r and A kept constant inside the integrals. This yields 

i 



a - 1 & _ A_» 
A °""2' E~^ T ° 



-1/2 



12) 



In ( V 1 + ^ 4 + r 2 ) - 2 (r 4 + l) 
In particular, Eqs. ( |T2"D predict that the DM soliton propagates steadily at anomalous 



average dispersion, /3 < 0, provided that Tq > (tq) ci ~ 0.301, at /5q = if Tq = (r, 2 



y'cr' 

2 ^ /V2\ /^2A 



and at normal average dispersion, /?o > 0, if Tq < (tq) ct . (To) min ~ 0.148 

In the case (3q > 0, Eq. fli~2"D predicts that the solution exists in a limited interval of 
the normal- PAD values, viz., 

< fa/E < » 0.0127. (13) 



Moreover, inside this interval Eq. (12) yields two different values of the minimum width 
Tq for a given (3qJ E (while in the anomalous- PAD region, r is a uniquely defined function 
of 1%/E) 0. On the basis of general stability criteria |16|, one can immediately conclude 



that the solution (i.e., DM soliton) corresponding to the larger value of tq is stable, while 
the one corresponding to the smaller tq is unstable. The border between the stable and 
unstable solitons corresponds to (3q/E = {/3o/E) mgx , and it is at Tq = (To) min ~ 0.148. 
Translating Tq into S according to Eq. (|^) (in particular, (tq ) min gives rise to 5 max ~ 9.75), 
we eventually conclude that VA based on Eqs. (|5|) and (|T0p predicts the following: 

• stable DM solitons at anomalous path-average dispersion if S < S CT ~ 4.79; 

• stable DM solitons at zero path-average dispersion if 5* = S CT ~ 4.79; 

• stable DM solitons at normal path-average dispersion if 4.79 < S < S ma _ x ~ 9.75; 



• no stable DM soliton if S > S max ~ 9.75 . 

Below, we will use, instead of Pq, the power normalized to that of the fundamental 
sech soliton having the same width as a given pulse ||, which is P = 4 ■ 1.12P (the factor 
1.12 is the ratio of the FWHM for the sech-shaped and Gaussian pulses). To further 
illustrate the properties of the solitons in the periodic-DM model, in Fig. la we show the 
normalized power P vs. the map strength S for different values of (3q, as predicted by 
Eqs. fljjD and (0) (the dependences are shown only in the region S < 9.75, where the 
solitons are expected to be stable). For comparison, Fig. lb shows the same dependence 
obtained from direct simulations of the full equation ([!]). The stars mark in Fig. 1 
particular solutions whose response to random variations of the fiber segment lengths 
will be displayed below. The curves in Fig. lb corresponding to normal PAD (/? > 0) 
terminate at points where the corresponding DM soliton becomes unstable. 

The comparison of Figs, la and lb shows that VA based on Eqs. (|6|) and fliTf ) yields 
quite acceptable results (for periodic DM) just in the range of small energies/powers, 
for which this approximation was devised. In particular, the VA-predicted S CI ~ 4.79 is 
different from but nevertheless close to the critical DM strength S CT ~ 4 which the direct 
simulations give for the small-power case. With the increase of power, the numerically 
found S CT grows, as is evident in Fig. lb. It is also noteworthy that the value S'max ~ 9.75, 
predicted by VA as the stability limit for the DM solitons, is indeed close to what is given 
by the direct simulations for small powers, see Fig. lb. 

At larger powers, there is a considerable discrepancy between VA and the direct numer- 
ical results. However, simulations displayed below show that, in the random-DM system, 
the soliton is strongly unstable at large energies anyway, so we are really interested only 
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in the small-energy range, for which the above VA is adequate. 



4 The randomly dispersion-managed system 



In the general case, when the pulse transmission is not steady (including the random-DM 
case), its evolution from a cell to a cell can be described in terms of a map, r — > r + 6r , 
A — > A + d~A . Because the changes are small, many iterations of the map, cor- 
responding to the propagation distance comprising many DM cells, may be approxi- 
mated by smoothed differential equations, dr /dz = 5t / +L^ + ^and dA /dz = 
5A / A straightforward calculation, using Eqs. (|TTD and (|i~0|) and taking 
into regard the normalization |/?_| = = 2 adopted above, leads to a final form of the 
smoothed equations: 



dr 

dz 



y/2E 



V -J^Tq , 1 

1 [L(-) + L«] 1 



+ 



1 



V 4 + 4Ag ' ^r 4 + 4[A + 2L(-) -2L(+)] 2 
1, 



V 4 + 4[A + 2L(-)] 2 ' 



(14) 



dA Q 

dz 



V2Et3 2A, 
[IH + £(+)] 1 

A + 2L< 



2 [A + 2L(-) - 
r 4 + 4Ag yV 4 + 4 [A + 2LH - 2L 



4 A + 2L^ i / r - 

L J =f - - In 2A + v/ro 4 + 4A, 

o 4 + 4[A + 2L(-)] 2 2 V 

^ In ^2 [A + 2LH - 2L« 

A + 2L ( - } 



+ \r ( 



- 4 + 4[A + 2L(-) -2L(+f) 



+ ln 2 



+ 



t 4 + 4[Ao + 2L(-)] 2 )}. 



(15) 



To check these equations, one can get back to the case of periodic DM, with ) 
Jji+i = 1/2, as per Eqs. @. In this case, a fixed point of Eqs. ([14]) and (|T5|) (dr /dz 



dAo/dz = 0) has exactly the same values of Ao and tq as given by Eqs. dl2|) . In the 
next section, we will display results of numerical integration of Eqs. (|T4"D and (|l^) for the 
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1L model of random DM defined above, in which the values = iS -* are picked up 
randomly from the interval (|5|). 

5 Numerical simulations of the variational equations 

Equations (f[4[) and (15) with random coefficients were numerically integrated, with initial 



conditions corresponding to the chirpless pulse, whose parameters were taken as per the 
fixed point ( |12"D of the allied periodic-DM model. The most essential single characteristic 
of the pulse propagation at given values of (3 and E is the rms cell-average pulse's width 
W, which we define as 

W = L~ l <j>W{z)dz (16) 

where the relation (§) and normalizations @) have been used. It is noteworthy that, in 
the case of periodic DM, the steady propagation regime corresponding to the fixed point 
(|T2| ) with A = —1/2 gives rise to the minimum rms width at Tg = l/v^3 (i.e., in the 
anomalous-PAD region, as l/\/3 > (tq) ct ~ 0.30). 

Simulations of Eqs. fll^D and fll5|) reveal that there are two different dynamical regimes. 
In the case when the soliton's energy is sufficiently low, i.e., one is indeed close to the quasi- 
linear approximation for which the above derivation is relevant, and PAD is anomalous 
or zero, i.e., (3q < (especially, if (3q = 0), the pulse performs random vibrations but 
remains, as a matter of fact, fairly stable over long propagation distances. In the case 
when the energy is higher, as well as when PAD is normal, (3$ > 0, the pulse demonstrates 
fast degradation and spreading out. 

To present the results in a more systematic form, in Fig. 2 we display the results 
for a typical case of the effectively stable propagation with small normalized power of 
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the soliton, P = 0.16 (in this case, E = 0.5), and anomalous PAD, (3q = +0.2 (which 



is 10% of the local dispersion). The simulations of Eqs. (R) and (|U|) were performed 
200 times, with the same initial conditions but different realizations of the random length 
set. Using numerical data for the 200 runs, we computed the evolution of (r ) and the 
averaged dependence (w(z)) ((...) stands for the averaging over 200 runs), along with 
the corresponding standard deviations. Both dependences are displayed in Fig. 2. 

It is clearly seen in Fig. 2 and in a number of similar plots not displayed here that, 
on top of the random vibrations, which are directly induced by random DM but are 
eliminated by averaging over 200 realizations, the soliton demonstrates slow (long-scale) 
dynamics. Systematic degradation of the oscillations and of the soliton itself takes place 
too, the degradation being slower for lower powers. For the particular case shown in 
Fig 2, the pulses remain certainly usable over the propagation distance ~ 100 DM cells. 



Further simulations of Eqs. (14) and (15) for propagation distances essentially larger 
than 1000 DM cells (not shown here; see an example in Ref. |J) show that, in fact, the 
sluggish spreading out of the soliton suddenly ends up with its blowup (complete decay 
into radiation). It is interesting to note that VA has predicted the same scenario for 



the evolution of the soliton in the early work |L4j for a model with regular sinusoidal 
modulation of the local dispersion: a long span of chaotic but nevertheless quasi-stable 
vibrations is suddenly changed by rapid irreversible decay. 

Fig. 3 shows a drastic difference in the soliton's evolution which takes place if the 
power is increased to P = 0.44, (the corresponding energy is E = 2.5), without changing 
parameters of the random-DM fiber link. In this case, which is typical for high powers, 
rapid decay of the soliton without long-scale vibrations is observed. Note that, since (tq) 
grows quite slowly, and the spectral width of the pulse is ~ Tq 1 , most of the pulse's 
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broadening in the temporal domain observed in Fig 3 is due to chirping of the pulse, 
rather than directly to a change in its spectral width. 

Taking f3 = 0, instead of anomalous PAD, radically improves the situation, as is 
seen in Fig. 4. For the high-power case with P = 0.44 (E = 3.6), which gave rise to 
rapid degradation of the soliton in the presence of the anomalous PAD, the pulse now 
survives over much longer distances, see Fig 4a. For lower values of the power, we observer 
still more robust propagation with zero PAD. Fig 4b shows the case P = 0.12 (E = 0.1), 
where the pulse shows slow dynamics induced by the random length variation in a strongly 
dispersion-managed link (S ~ 4.8) with very little degradation over 1000 DM cells. 

The PAD's value (3q = turns out to be a point of a sharp optimum: taking any 
tangible small normal value of PAD, /5 > 0, we always observed rapid decay (without 
long-scale oscillations) of the soliton at virtually all the values of the energy, see for 
instance Fig 5, which displays the case with P = 0.27 (E = 2.5) and /3 = 0.02 (1% of 
the local dispersion). In this case the soliton character of the pulse is lost quickly. The 
broadening of the pulse in the spectral domain is enhanced by its spectral broadening, 
manifested in the decrease of To- The broadening rate, however, is slower than in Fig 3a 
due to the lower magnitude of PAD. 

6 Direct simulations 

It is necessary to compare the predictions of VA with direct simulations of Eq. [I]. The 
pattern of the simulations was the same as in the previous section, i.e., simulations were 
performed for 200 different realizations of the random length set, in order to evaluate 
averaged evolution of the pulse's parameters. First, the shape of the DM soliton in the 
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allied periodic-DM link was numerically determined, to be used as the initial configuration. 
The randomness was the same as above, i.e., the lengths were uniformly distributed 
between 0.2 and 1.8 of the average length. The general trends predicted by the VA are 
confirmed by the direct simulations: lower energy and anomalous or, especially, zero PAD 
enhance the stability, see details below. 

In direct simulations, the width must be defined with special care. The usual FWHM 
definition may be misleading in the case of random DM, as the soliton can sometimes split 
(see below); besides that, this definition ignores the presence of a radiation component. 
Therefore, we adopted an integral definition, with which the width is a size of the temporal 
region on both sides of the soliton's center that contains 76% of the net energy. For a 
pulse with the Gaussian shape the width defined this way coincides with FWHM. We 
will assume that optimum dispersion compensation can be applied at the receiver, i.e., 
any linear chirp is removed from the pulse by a dispersion compensating element installed 
before the receiver. Note that in the framework of the above VA, r represented the 
pulse's width at the chirp-free points. 

Starting with anomalous PAD, in Fig 6 we show results for the same case for which 
results obtained by means of VA were displayed in Fig 2, i.e., the magnitude of the 
PAD is 10% of the local dispersion and the power is low, P = 0.18. This time, we 
show (by bold curves) not only the averaged evolution, but also (by thin curves) the set 
of the particular evolutions corresponding to different realizations of the random length 
set. Comparing Figs. 6 and 2, we conclude that the averaged results are qualitatively 
similar. In particular, internal vibrations are present in the first 100 cells in both cases, 
and the systematic temporal-domain broadening of the pulse takes place in both cases, 
although the broadening rate is overestimated by VA (hence, the full numerical results, 
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predicting slower degradation of the soliton, seem considerably better for the applications 
than the less accurate results generated by VA). The latter discrepancy is, most plausibly, 
accounted for by radiative losses that are ignored by the VA. 

Since we allow for dispersion compensation at the receiving end, this can change the 
link's PAD. The effective PAD, corrected with regard to the optimal receiving-end disper- 
sion compensator, may therefore be regarded as a function of the propagation distance. 
The inset to Fig. 6 shows that the thus redefined effective PAD stabilizes at the value 0.2, 
which turns out to be the same as for the unperturbed DM soliton. This means that no 
post-transmission dispersion compensation is needed in order to receive a chirp-free pulse; 
the dispersion of the link is balanced by its nonlinearity. For higher powers, however, this 
is no longer the case, see an example for P = 0.47 in Fig. 7. The pulse broadens and 
decays fast and the optimum end dispersion drops rather than stabilizing, i.e., the pulse 
acquires a chirp during the propagation and extra post-transmission dispersion compen- 
sation will be necessary. Compared to the variational results, displayed for the same case 
in Fig 3, we conclude that the direct simulations of Eq. ([!]) reveal a rapid decay of the 
pulse shape. Moreover, detailed e4xamination of the pulse shapes generated by the direct 
simulations shows they often develop a multi-peak structure, so that it can sometimes be 
hard to identify the pulse proper. This is why Fig. 7 shows the evolution for the first 50 
DM cells only. 

It is necessary to mention that essentially longer direct simulations of Eq. ([I]) (not 
shown here) demonstrate that the gradual decay of the soliton is, at a final stage, suddenly 
interrupted by splitting of the residual pulse into two smaller ones. Recall that, although 
VA by definition cannot predict the eventual sudden splitting of the pulse, it does predict, 
as it was mentioned above, something similar, viz., sudden fast decay of the soliton after 
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a long stage of chaotic vibrations. It is noteworthy that gradual evolution of the soliton in 
the model with sinusoidal periodic modulation of the dispersion also ends up with sudden 
splitting, which approximately corresponds to the sudden decay predicted by VA in the 



same model |TR [19]. 



To further illustrate the internal dynamics of the pulse, we picked a case where PAD 
is anomalous and equal 1% of the local dispersion, and the power is P = 0.1. Fig. 
8 shows the average trajectory on the dynamical plane, where the coordinates are the 
width and chirp of the soliton at the midpoint of the anomalous-fiber section. The chirp 
is in this case defined as the second derivative of the phase of the pulse taken at its 
center. The trajectory demonstrates a quasi-circular motion around a center which is 
drifting to the right, towards a broader pulse. The trajectory clearly demonstrates the 
nonlinear character of the system, cf. the inset to Fig. 8 which shows the trajectory 
if the nonlinearity are dropped. Note that a stationary periodic-DM soliton would be 
represented by a single point in this plane, and a perturbed (nonstationary) soliton in 
the periodic-DM system would trace a circle around a fixed center shifted to a broader 
pulse, as compared to the stationary DM soliton [pOf] . The random length variation is 
a permanently acting perturbation which generates the persistent drift of the center in 
Fig. 8. 

Numerical simulations at zero PAD also agree well with the results generated by VA: 
zero PAD provides the slowest pulse degradation in the random DM link. Fig. 9 shows 
the evolution of the pulse width when /3 = and the power is P = 0.15. Almost no 
broadening (< 5%) can be observed, on the average, after having passed 200 DM cells. 
The broadening does increase with the power, but remains, as predicted by VA, much 
less than in the case of nonzero PAD. For instance, taking P = 0.45 results only in a 
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doubling of the width on the propagation distance equal 200 DM cells. Comparison to 
the variational results for the same case displayed in Fig. 4 again shows that full numerical 
results are qualitatively similar but somewhat better, predicting slower degradation of the 
pulse. One of the reasons for this may be that the radiation remains trapped within the 
pulse. 

Lastly, for the case of normal PAD the direct simulations confirm the prediction of VA 
according to which the solitons are least stable in this case. This is partly explained by the 
fact that there is, effectively, a minimum (threshold) power necessary for the existence 
of stable solitons at normal PAD even with periodic DM, see Fig lb, and the larger 
power always stimulates the degradation of the pulse. Detailed examination of the direct 
numerical results for the normal-PAD case also shows that the pulse's spreading out is 
boosted by its spectral broadening during the initial stage of propagation. As a result, 
even at very weak normal PAD, /3 = 0.02 (1 % of the magnitude of the local dispersion), 
the pulse stays intact for no more than 10 DM cells. 

Conclusion 

In this work, we have put forward a model of a long optical communication line subject to 
random dispersion management (DM). The line consists of alternating fiber pieces with 
anomalous and normal dispersion, whose lengths are picked up randomly from a certain 
interval, while the absolute values of the dispersion coefficients in both pieces are always 
equal. By means of the variational approximation, we calculated small changes of param- 
eters of the propagating quasi- Gaussian pulse per one DM cell, and then approximated 
the evolution of the pulse passing many cells by smoothed ordinary differential equations 
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with randomly varying coefficients. The equations were solved numerically, and results 
were presented as the dependences of the pulse's mean width, and average deviation of 
the width from the mean value, on the propagation distance. Averaging the results over 
200 different realizations of the random length set removes rapid oscillations of the width 
and reveals slow long-scale dynamics of the pulse, frequently in the form of long-period 
oscillations. The second part of the work is based on direct numerical simulation of the 
same model. Comparing the results, we have concluded that essential features of the 
soliton dynamics are the same in the variational approximation and in direct numerical 
simulations: the propagating soliton is most stable in the case of zero path-average dis- 
persion (PAD), less stable in the case of anomalous PAD, and least stable in the case 
of normal PAD. The soliton's stability also strongly depends on its energy, so that the 
soliton with small energy is much more robust than its large-energy counterpart. 
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FIGURE CAPTIONS 

Fig. 1. The normalized power vs. the map strength for stationary solitons in the periodic 
DM system: (a) the analytical result, Eq. (|T2"D, produced by the variational approxima- 
tion; (b) direct numerical solution of the NLS equation ([]]). Each line corresponds to a 
constant value of PAD: from left to right, (3q = —0.2, —0.02, and 0.02. The stars mark 
particular cases for which response to random length variations are further displayed in 
Figs. 2 through 9. 

Fig. 2. The cell-average pulse width (top) and minimum-width parameter r (bottom) 
vs. the propagation distance, generated by numerical integration of the variational equa- 
tions (|T4| ) and ([15]) for the power P = 0.16 and (3q = —0.2 (anomalous PAD). The mean 



values (solid curve) and standard deviations from them (dashed curves) are produced by 
averaging over 200 different realizations of the random length set. 
Fig. 3. The same as in Fig. 2 but for higher power, P = 0.44. 

Fig. 4. Evolution of the cell-average pulse width for zero PAD propagating over 1000 
DM cells, in the cases of high power P = 0.47 (top) and low power P = 0.1 (bottom). 
The mean values (solid curve) and standard deviations from them (dashed curves) are 
produced by averaging over 200 different realizations of the random length set. 

Fig. 5. The same as Fig. 2 but with normal PAD, f3 = 0.02, and P = 0.27. 

Fig. 6. The pulse width vs. the propagation distance for anomalous PAD, /3q = —0.2, 
and low power, P = 0.18, generated by direct simulations of Eq. (|1|). Shown are both the 
mean values (solid curve) and standard deviations from them (dashed curves) as produced 
by averaging over 200 different realizations of the random length set, and grey curves 
corresponding to particular realizations of the random set. The inset shows the effective 
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PAD corrected with regard to the optimal dispersion compensation at the receiving edge 
(average of 200 simulations). 

Fig. 7. The same as in Fig 6. but for high power, P = 0.47. 

Fig. 8. Results of direct simulations of the NLS equation ([I]) with anomalous PAD, 
(3q = —0.02. The plot shows the pulse dynamics in the dynamical plane whose coordinates 
are the chirp at the center of the pulse and its width, evaluated at the midpoint of the each 
anomalous-fiber section. The inset shows the same but in the absence of the nonlinear 
term in Eq. (1). 

Fig. 9. The same as in Fig. 6 for the case of zero PAD and P = 0.15. 
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